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Interfacial structure and correlation functions near a two-dimensional (2D) wedge filling transition 
are studied using effective interfacial Hamiltonian models. An exact solution for short range bind- 
ing potentials and results for Kratzer binding potentials show that sufficiently close to the filling 
transition a new length scale emerges and controls the decay of the interfacial profile relative to 
the substrate and the correlations between interfacial positions above different positions. This new 
length scale is much larger than the intrinsic interfacial correlation length, and it is related geomet- 
rically to the average value of the interfacial position above the wedge midpoint. The interfacial 
behavior is consistent with a breather mode fiuctuation picture, which is shown to emerge from an 
exact decimation functional renormalization group scheme that keeps the geometry invariant. 
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I. INTRODUCTION 

Fluid adsorption in wedge and cone-shaped non-planar 
geometries has attracted much attention in the last few 
years P, Q, S S E| ■ Geometry plays an important role 
in the surface phase diagram, and new phase transitions 
as the filling transition arise. Thermodynamic considera- 
tions 0,|7|,|3] predict that the gas- liquid interface unbinds 
from the wedge before the wetting temperature corre- 
sponding to the substrates. So, the wedge is completely 
filled by liquid for temperatures higher than the filling 
temperature Tf < T^, where Tf is given by the condi- 
tion: 



OiTf) 



(1) 



and 9{T) is the temperature-dependent contact angle of 
a liquid drop on the planar substrate. Capillary wave 
models show that the filling transition can be critical 
even though the wetting transition corresponding to the 
substrate is first order, and that interfacial fluctuations 
are enhanced with respect to the wetting case 0, Q ■ For 
the 2D wedge filling transition in shallow wedges charac- 
terized by an small angle a respect to the x axis (see 
below), there exists a remarkable covariance relation- 
ship between the wedge midpoint probability distribution 
function P^{lo) in the filling fluctuation regime and the 
planar f— point probability distribution function P^{lo) 
characteristic of a strong-fluctuation regime critical wet- 
ting transition: 



(2) 



where 9 is the contact angle of the liquid droplet on the 
substrate. This expression establishes a connection be- 
tween two apparently unrelated phenomena, the deep ori- 
gin of which is still elusive. The covariance relationship 
has been observed also in acute wedges 9], Ising model 
exact calculations lOj and computer simulations f L] . Al- 
though the covariance relationship is restricted to the in- 
terfacial behavior above the wedge midpoint, some other 



quantities like the local susceptibility, which is related to 
the 2— point correlation function, also showed a modified 
covariance relationship . Consequently, it is interesting 
to see if the covariance extends to higher-order probabil- 
ity distribution functions. 

In this paper we study the structure of the interfacial 
profile and correlations for 2D wedge filling phenomena. 
Exact results for the capillary wave effective Hamiltonian 
theory in the filling fluctuation regime are obtained as an 
extension of the analysis presented in Ref. j^^- The ex- 
act results show the appearance of a new length scale 
across the wedge close to the critical filling transition. 
This scale controls the decay of the interfacial profile, 
local roughness and correlations, and is related geomet- 
rically to the wedge midpoint average interface position. 
For the local properties, we found a very interesting re- 
lationship between the wedge 1— point probability distri- 
bution function and the corresponding functions in the 
planar geometry, which can enlighten the origin of the 
wedge covariance. 

Regarding the two-point correlation functions, we 
found a confirmation in the scaling limit of the breather 
mode picture ^ A\ , which states that the interface is ef- 
fectively infinitely stiff in the filled region and is driven 
by fluctuations of the wedge midpoint interfacial posi- 
tion, i.e. critical effects at 2D wedge filling arise simply 
from local translations in the height of the flat, filled in- 
terfacial region. 

Finally, we explain the critical behavior of the fill- 
ing transition in the functional renormalization group 
approach. As the geometry is fundamental in the un- 
derstanding of the critical filling transition, we choose 
a scheme that leaves the wedge geometry invariant. 
We show that the breather mode picture emerges as a 
straightforward consequence. The predictions for the 
critical behavior are in complete agreement with exact 
solutions. 

Our paper is organized as follows. In Section ^] we 
describe the continuous transfer matrix formalism and 
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FIG. 1: Schematic illustration of a typical interfacial configu- 
ration in the wedge geometry. The relevant correlation length 
scales and ^±{x) are also highlighted. Other notation is de- 
fined in the text. 



the definition of the wedge rt— point interfacial probabil- 
ity distribution functions. We apply this formalism to 
the case of contact binding potentials in Section ITTll and 
in particular calculate analytically the 1— point proba- 
bility distribution function and the 2— point correlation 
functions. Some results for Kratzer binding potentials 
will be presented in Section IWI In Section we analyse 
the breather mode picture and derive a relation between 
two important scaling functions. Section IVII is devoted 
to the development of a renormalization group theory of 
2D critical filling transition, which requires a generaliza- 
tion of previous approaches for critical wetting. A brief 
conclusion is presented in Section IVIII 



II. THE FORMALISM 

Consider a two-dimensional wedge formed by the in- 
tersection of two equal planar substrates at angles zta to 
the horizontal (see FigQ]). We suppose that the wedge is 
in contact with a bulk vapor phase at saturation condi- 
tions, i.e. in equilibrium with the liquid phase, and the 
substrates preferentially adsorb the liquid phase. Our 
starting point is the effective interfacial Hamiltonian for 
shallow wedges: 



pH[l] 



X/2 



-X/2 




Wiyix)~a\x\)\ (3) 



where y{x) is the interfacial local height respect to the 
horizontal, X is the interfacial horizontal length, ksTYi 
is the interfacial stiffness, kBTW{l) is the local binding 
potential and f3 = l/ksT. We impose periodic boundary 
conditions at the ends, i.e. y{—X/2) — y{X/2). 

Defining the local relative height between the vapor- 
liquid interface and the substrate l{x) = y{x) — a\x\, 



Eq.Q can be rewritten as 0: 
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\dx J 



dl 



Eal ^ ) {2H{x)-l) + W{l{x)) 



(4) 



where H{x) is the Heaviside step function. Integrating 
by parts to eliminate the term proportional to (dl/dx), 
the effective Hamiltonian can be expressed as 



(3H[l] = X 



x/2 

-X/2 



2T.al{X/2) - 2Y.al{0) 



(5) 



The first two terms in the equation are irrelevant con- 
stants for the interfacial properties in the wedge, the 
third one is the origin of the boost factor that decreases 
the pinning effect of the binding potential 0, and the 
fourth one corresponds to the effective Hamiltonian of an 
equivalent planar interface problem. As the probability 
distribution of an interfacial configuration is proportional 
to exp{—f3H) we can relate the wedge and planar prob- 
ability distributions in a straightforward way. In partic- 
ular, the n— point wedge correlation functions can be re- 
lated to (n + 1)— point correlation functions in the planar 
case by adding the wedge midpoint position. However, 
the presence of the boost factor will alter significally the 
behavior of the wedge correlation functions with respect 
to their planar counterparts. 

Our approach is based on a standard application of 
transfer matrix methods |13|. The partition function 
ZTr{li,h,xi,X2) of the interface with fixed endpoints 
(xi,Zi) and (x2,?2) with X2 > xi in presence of a pla- 
nar substrate is defined as the following path integral: 



Zt,{IiJ2,Xi,X2) = ZT,ili,l2;X2 - Xi) = 

S /dT ' 
2\dx. 



= Vlcxpl - 



dx 



W{1) 



(6) 



The partition function Eq. ^ is the solution of the 
following Schrodinger equation: 
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with the initial condition: 



Z^{h,l2;x) =0 (7) 



z^{h,l2;0) = Sil2-h) 



(8) 



where d{x) is the Dirac delta function. Formally, the 
partition function can be expressed as: 

Z^{h,l2■,x) = J2^Uh)4'^{l2)eM-E^x) (9) 
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where ipH) and Ei are the eigenfunctions and eigenvalues 
of the time-independent Schrodinger equation: 



1 

"2E 



(10) 



with appropriate boundary conditions. In the thermo- 
dynamic Umit ^ exp{—(3fX) as X oo, where 
/?/ = E(cos 9 — 1) is the excess free energy per interfacial 
length. Consequently, Eq. Q implies that Eq = /3/, so 
that in the low contact angle limit, Eq « — E0^/2. 

The n-point distribution functions can be obtained in 
terms of Zt,[Ii, I2] x) as: 



: n) = lim 



ni=0 ^Tri^'h h+l', Xi^i — Xi) 



i=l 

where i = [l^Xi), Xn+i = -xq = X/2, and /q = In+i = 
lx/2- For n = 1, P^{i) = iV'ol^OP- From Eqs. and 
@ it is clear that if the distance between two subsets 
{xi, . . . , Xm\ and {xm+i, . . . , a;„} is much greater than 
the planar correlation length = l/(£'i — E^) (with Ei 
the first excited state eigenvalue), the distribution func- 
tion factorizes and the two subsets become uncorrelated. 

The n-point wedge distribution functions /\o(l; . . . ; n) 
can be expressed, in general, in terms of {n + 1)— point 
planar distribution functions. So, for a set {x_.,„ < 
x_i < < < . . . Xn}-, they can be expressed as: 

J2T,alo 



< 



Pni{-m; . . . ;n) 
-l;0;l;...;n) 

where (n|/(/)|m) 
. . . < 
pier: 

P^(l; 



din 



|0) 

P^(-l;l)P^(-m;. 



(12) 



< Xn, the expression of ^^,(1 
n) 



j;^diMi)f{i)rm{i)- If < X, < 

; n) is slightly sim- 
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A similar expression is found if < ... < Xn < 0. 
Finally, if a; = is included in the x set, the wedge n- 
point distribution function reduces to: 



Pwi-m; . . . ;n) 



-P,(-m;...;n) (14) 



Although this approach is general for arbitrary binding 
potentials, we will restrict ourselves to some special cases. 
The first case will be contact potentials, in which W{1) — 
for Z > 0, W{1) = +00 for Z < and at the wall the 
eigenfunctions fulfill the boundary condition [l^ : 



d_ 
di 



\TLi>{l) 



(15) 



where r is proportional to the deviation from the crit- 
ical wetting temperature. For t > the contact an- 
gle is related to r via t — HO 13]. These potentials 
can be understood as the limiting case of a square-well 
binding potential when the well width tends to zero. Its 
importance is threefold. First, this case corresponds to 
the filling fluctuation regime, that previous studies show 
to be the relevant one for potentials which decay faster 
than \/l. Secondly, there is an analytical expression for 
given by: 



Z^{li,l2;x) ^ \l — I e 2x +e 2, 



.ew-^(..^.)erfc(,/ — (/,+y-rJ^ (16) 



Finally, this case can be compared to more microscopic 
results, like the exact solutions of the interfacial proper- 
ties of the corner filling of an Ising model. 

Another interesting case is the Kratzer binding poten- 
tial [l3: 



(17) 



where (j) = {1 + \/l + 8Ew)/2 and we assume Dirichlet 
boundary conditions at the origin. Previous studies indi- 
cate that this class of binding potentials corresponds to 
the marginal case between the mean-field and fluctuation 
dominated regimes for the critical filling transition. The 
Laplace transform of Z-„[li,l2]x), Zt;{Ii,12, E), is given 
by Hi: 



Z^{h,l2,E) = / dx e''^Z^{h,l2;x) 
Jo 



fr 



r [20] 



w 



(18) 



•-0 



where Eq = —I]6'^/2, /> = max(/i,Z2), l< = niin(Zi,^2), 
r(a;) is the gamma function, and finally Mi^^m{z) and 
W^./miz) are Whittaker functions, related to confluent 
hypergeometric hmctions. 



III. EXACT RESULTS FOR CONTACT 
BINDING POTENTIALS 

In this Section we will obtain and analyze some rel- 
evant wedge distribution functions for contact binding 
potentials. In particular, we will revisit the 1— point dis- 
tribution function (considered previously by our group 
[T^l and the 2— point height-height correlation function 
between the midpoint and any other interfacial posi- 
tions. Related quantities as the average interfacial profile 
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{l{x))w, the local roughness ^±{x) and the correlation 
length across the wedge £,x (see Fignj will be also ob- 
tained. 

Some results are already known for the 1— point distri- 
bution functions. The probability distribution function 
for the midpoint a; = interfacial height is given by 0|: 



P^Go;0,a)^P»(/o,O) = 2E(0-a)e-2S(«-")'« (19) + 



that verifies the remarkable covariance relationship Eq. 
©■ 

For arbitrary a; > the 1— point distribution function 
has the expression^Q : 



V {l'^{x))w - {K^))^^ where {l'^{x))nj is given by: 
1 



1 



a) S0 




^^erfciy:j:(.-2a) 



1 



6l2 6l2 \ 9x a 



— a \ X 



2fl2 




(23) 



Ex „ / E 



P^(/,a;) = Eee-^^^^erfcJ^-^— (?+-^— / 
+ ^{9 - a)e2S(a-e)/e2S-(--e) 




- Eae-2^"'e2^"-("-«) 

X erfc(y^(^-2a) + y^z) (20) 



For general n, the following expression can be obtained 
by induction: 



f_ \ 2'Sax{a-e) 




+Qn{x)eTk I \24) 



For X < 0, we have the symmetry Pyj{l, x) — Pw{l, ~x), 
so hereafter we will consider only the case a; > 0. 

The moments {l^(x))w can be obtained after some al- 
gebra. The average interfacial position profile reads: 



1 

2E^ 



X st)^ 

-e 2 



27rE 



a a 



4Eg ^^^MV~^ 



2a) 
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4E^ 
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/E6I2 



-x (21) 



The wedge excess adsorption Fu, measured with respect 
to the planar case can be obtained as: 



F,„ = 2{pi - pg) (^{l{x))^ - dx 
pi- Pg \ 1 1 



2E2 



61(61 - a)2 



613 



(22) 



where (/")7r = n!/(2E6')" and P„(x) and Qn{x) are poly- 
nomials in X of order n — 1 and n, respectively. 

These expressions are only valid if > a (for smaller 
values of 9 the interface is unbound from the wedge). For 
X ^ 0, Eq. if^ reduces to Eq. ifT^ll . On the other hand, 
for |a;| oo, Pw{l,x) decay to P^{1) = 2E6' exp(-2E6'/). 
However, the scale over which this decay occurs depends 
on the value of a. If 6* > 2q!, this scale is the planar cor- 
relation length = 2/E6'^. However, if a < 6* < 2a, the 
decay length is = l/2Ea(0 — a) (our notation differs 
slightly from the one used in Ref. [131 )■ Note that S,f is 
always larger than ^|| , and diverges on approaching the 
filling transition. On the other hand, £,f is related ge- 
ometrically with the wedge midpoint average interfacial 
height via S^p — {l{^))w/oi ~ (Z(0))^/tana for small a. 

It is amusing to note that Eq. H2U|) verifies the following 
differential relation: 

PUi,^) + ^F (^^4^)=LAi,^) 



dx 

oo 



1 d 

= PAl) + -^Q^ J dlokPAlo, 



0;l,x) 



(25) 



where L.,j[l,x) is for contact binding potentials: 



where pg and pi are the coexistence densities of the va- 
por and liquid phases, respectively. Close to the filling 
transition {9 a), F^, 2{pi - pg){l{Q))'^/a. 

The roughness profile ^±{x) (see Fignj is defined as 



L^il,x) = J:9e 



-2S«'erfc(-A/— 6 



2x , 



^'^'i-x'^ 



(26) 



Note that the RHS of Eq. H25|l depends only on the pla- 
nar properties and, consequently, is independent of a. It 
can be shown that Eq. (|25|l is obtained for any binding 
potential if the LHS is expanded in powers of a and trun- 
cated at the lowest order term, which is independent of a. 
Consequently, this differential field equation implies an 
infinite hierarchy of integro-differential relationships for 
the 2— point planar correlation function. Alternatively 
Eq. H25|l provides an elegant route to the calculation oi 
any moment of the interfacial height. Multiplying Eq. 
(|25|l by arbitrary power of / and integrating over all the 
possible values of I, the following differential equations 
are obtained: 



d(r(x)). 



dx 



dl Z"L^(Z,x) 





1 d 



in^ + Q^iiion^))^ (27) 



where {. . ■)uj and (. . .)7r mean the average with the wedge 
and planar distribution function, respectively. The RHS 
of Eq. H27(l depends only in the planar distribution func- 
tions, and consequently decays to {1"^)^ for distances 
larger than ^y. Close to the fihing transition, S^p ^ ^||, 
and we can approximate Eq. (|27|l for x > by: 



dx 



(28) 



that has as a solution (Z"(x))„ « (r)^ + ((/"(O))^ - 
(?"'),r) exp(— Taking into account that (?"(0))u, » 
(P),r close to the filling transition, the approximate 
solution can be simphfied even further to (/"(x))^, « 
(^"(0))u, exp(— x/^f) (which is equivalent to set {1^)tt = 
in Eq. (|28|l ). These findings are obviously in agreement 
with Eq. H24|l and the asymptotic behavior of Pw{l,x) 
for large x and 6 <2a 12]. 

It is interesting to note that the moments obtained 
from the actual 1— point distribution function are the 
only solutions of Eq. (|27|l that: (a) decay exponentially 
within a length scale lor < a/9 <^ 1 and a; ^ oo; 
(b) are analytical as a function of a for < a < in 
particular at the disorder point. The existence of the 
relationship Eq. H25() , from which covariance for the mo- 
ments of the interfacial position profile at a; = can be 
inferred provided the (a) and (b) regularity conditions 
are fulfilled, leads us to speculate on the existence of a 
hidden symmetry of the hamiltonian that explains wedge 
covariance. However, the nature of such a symmetry (if 
any) is completely unknown. 

In the mean-field approximation, the average interfa- 
cial position profile for binding potentials characterized 
by a critical exponent ~ fulfills the following gener- 
alized covariance relationship |T5| : 



l{x) = 



dl{x) 



dx 



(29) 



where l{x) represents the (averaged) interfacial position 
at cc, and It:{9) is the planar (averaged) interfacial po- 




FIG. 2: Comparison between (l{x))w obtained in Ref. pHl 
by Ising model computer simulations for boundary magnetic 
fields h/ J = 0.595 (circles), h/ J — 0.597 (squares) and h/ J — 
0.599 (diamonds); and the approximation given by Eq. igSJ 
(continuous lines). The Ising model parameters are: a = 
7r/4, the temperature T = Tc/2 and the bulk magnetic field 
Hbuik = 0. The boundary magnetic field at the critical filling 
is he/ J — 0.606. The lengths \x\ and {l{x))jn are measured in 
lattice spacing units. See text for explanation. 



sition for a given contact angle 9. Making the substi- 
tution l{x) {l{x))u], it is clear from Eq. that 
this extended covariance is not verified for x (even 
asymptotically when a; — s- or — > oo). However, it is 
remarkable that there exists an analogous to Eq. H29|l . 
given by Eq. (|27|l for n = 1. 

To finish our discussion about the 1— point distribution 
functions, we compare our results with computer simula- 
tions of the 2D Ising model nj . Close to the filling tran- 
sition point, we expect that the approximate solution to 
Eq. (|28|l for n = 1 will be generalized for arbitrary a to: 



(0^ 



COS a 



(30) 



where now £^p is defined as (Z(0))u,/ tan a. We have tested 
this approximation with the simulation results reported 
in Ref. 11] (see Fig. |2J). The symbols correspond to the 
simulation data obtained for an square 64 x 64 Ising lat- 
tice with zero bulk magnetic field and boundary magnetic 
fields +h for the boundary rows ending at the lower left 
corner, and —h for the remaining boundary rows. In this 
geometry, a = n/A. The temperature is set to T = Tc/2, 
where Tc is the bulk critical temperature. For this tem- 
perature and a the critical filling transition occurs at 
he/ J = 0.606. Fig. 121 shows the computer simulation 
results for h/J = 0.595,0.597 and 0.599. We have no 
direct estimation of {1)tt. However, we have obtained 
(Z)^ by fitting the simulation data with \x\ < 16 lattice 
spacings (in order to minimize the effect of the upper 
left and lower right heterogeneous wedges) to the Eq. 
(|^ . The best fitting values are, in lattice spacing units, 
(Ott = 0.314,0.335,0.436 for h/J = 0.595,0.597,0.599, 
respectively. As it can be seen, the fitting to the Simula- 
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tion data is quite good, despite the crude approximations 
involved in Eq. (|3U|) . 

Now we want to characterize the 2— point correlations, 
in particular the correlations between the interfacial posi- 
tion above the wedge midpoint and the corresponding to 
an arbitrary x, which are given by the following function: 



2S I da 



{iix)m)^-(iix))^(i{o)). 



(31) 



Substituting Eq. ^ into Eq. we obtain: 

{iix)m) 



1 / fl2 



X (2a 



27rE 2Ea(6'-a) 

2 ■ 




{6- 2a) 132) 



This function decays exponentially to zero for large x. 
However, the characteristic correlation length (see Fig. 

depends on a: it is f || for 9 > 2a and if a < < 2a. 
Consequently, the disorder point not only introduce a 
new length scale for the average interfacial profile, but 
also for the interfacial fluctuations. 



IV. RESULTS FOR THE KRATZER BINDING 
POTENTIALS 

The Kratzer binding potential (see Eq. (|17|l 'l is the 
marginal case between the filling mean-field and filling 
fluctuation regimes. For such potentials the wedge mid- 
point probability distribution function also obeys wedge 
covariance Eq. (0): 



mo - «)] 



T[2(f>+1] 



exp(-2E(6l - a)lo) 
= P^{lo;0-al33) 



It is possible to extend the transfer analysis and obtain 
exact results for other quantities of interest. Consider, 
for example, the 1— point probability distribution func- 
tion Pm{l,x). The Laplace transform Pw{l\E) can be 
expressed as: 



'0 r[20 + i] 

X exp{-j:ei)Z^{lo,l,E-T,e^/2) (34) 



where ZttHq, I, E) is given by Eq. (|I8|I . This reduces to: 



P^{l-E) 




eV[2(t) + l]T[2(t)\ 



0gE(2a-e)/ow 



kT [(/)(1 - k)] 



/2Y.9lo 



^2I]6'/c 



70 S(2a-6/);o 



M. 



f2Y.ei\ 



(2Y.ei\ , ^ (2Y.eh\ 



(35) 



where k = l/^/l - 2E/Y.e'^. The poles of Pw{l]E) in 
the E real positive semi-axis are the characteristic in- 
verse length scales across the wedge of Pw{l,x). Since 
the second integral is over a finite interval and the in- 
tegrand does not diverges in that range, no new length 
scale emerges from it. For the first integral, we take into 
account that [T^ : 



x'^ exp{—px)WK.^f^{ax)dx 



r[ii + i^ + l/2]r[j^ - /i + I/2]a^+3 



r[i/-K+l](p + a/2)''+'^+3 

1 1 

— , u — K H 

2'^ 2 



X2F,[^, + ly+l,^l-K+l■,ly- (36) 



P+2 

where 2^1(0, b, c; x) is a hypergeometric function, li 9 > 
2a, the integral does not introduce any new characteristic 
length . However, for a < 6 < 2a a new singularity 
emerges for E(6l-2a)-hS6l/K = 0, i.e. E = 2Sa(6l-a) = 
1/S.F- Remarkably, has the same expression as for 
contact binding potentials, and is proportional (but not 
equal) to {l{0))w/a. 

From this it follows that the non-thermodynamic sin- 
gularity occuring at = 2a mentioned in the previous 
Section is not specific to contact potentials. A simple ge- 
ometrical argument given in Ref. explains why. The 
most relevant interfacial fluctuations are those where the 
interface leaves the substrate with a contact angle 6 (rel- 
ative to the tilted wall) at an arbitrary substrate point. 
If > 2a, the other side of the wedge does not play any 
role and we can anticipate that the only length scale that 
controls the 1— point distribution decay is ^|| . However, 
if 6* < 2a, the interface will eventually reach the other 
substrate, and consequently we can expect the geometry 
to play an important role leading to the emergence of 
a new length scale. Formally, this non-thermodynamic 
singularity occurs when the following integrals that arise 
from the spectral expansion of Zj^{li, l2]x): 



V'o(0exp(2Sa/)V';(0 



(37) 
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become ill-defined. There, tpp{l) are the scattering eigen- 
states with eigenvalues E = p^/2Y, and V'o(0 is the 
ground eigenstate. A straightforward WKB asymptotic 
analysis for the eigenfunctions shows that, for p ^ 0, the 
integrals given by Eq. H37I) become ill-defined for <2a 
for quite arbitrary choices of binding potential. 

As 6 /a decreases, exceeds the intrinsic interfacial 
length scales i/{Ei — Eq), and becomes the true corre- 
lation length across the wedge £,x at an another disor- 
der point when = ^|| (recall that = for 6 /a 
larger than the value at the disorder point). For the case 
of contact binding potentials both non-thermodynamic 
singularities occur the same value 9 = 2a. However, 
in general, the non-thermodynamic singularities are dis- 
tinct provided there are at least two bounded eigenstates 
of Eq. (|10|) . For the pure Coulomb case (0 = 1) the 
second disorder point occurs at 6 = Aa/3. 

Close to the new singularity $,p^ we found that: 



P^{l-E) 



1 



S6 



(38) 



so Pw{l,x) behaves asymptotically for large values of x 
as x^"-" exp(— x/^i?), provided that ^|| < S,f- 

A field equation analogous to Eq. (^5)1 can be found 
for Kratzer potentials. Transfer matrix calculations for 
arbitrary binding potentials lead to the relation: 



(0|e^^"'|0) 



e-a 



f dP^{l,x) 
\ dx 



dlo ^ty?,^ P^{Iq,0;1,x)\ = LAI, x) 



^OMk) 



' dlo j'ff, M io,0;l,x) (39) 



where ^/'o(^o) = V'oCo) exp(S6'/o) and iP'q{Io) is its deriva- 
tive respect to Iq (recall that "00 Co) is the ground state 
eigenfunction) . For Kratzer potentials, V'o(^o) Iq, so 
Eq. (|39l) can be expressed as: 



dPw{l,x) 



d_ 

da 



a J da 

dPw{l,x) 
dx 



(40) 



As for the contact binding potential case, some interest- 
ing quantities can be evaluated from this expression. For 
example, the wedge adsorption is found to be: 



r^-(20 + i)((/. + i)r' 



CP 



(41) 



where r*^^ is the adsorption corresponding to the contact 
binding potential Eq. H22|) . 



V. THE BREATHER MODE PICTURE 

In order to understand the origin of the new correlation 
length we identified in previous Sections, we recall 
the definition of the 2— point distribution function for 
X2 > xi > 0, Eq. (fT!^ . This expression can be written 
in the following way: 

P^{h,X2\li,xi) = PAl2,X2\li,xi) = PAl2,X2 - a;i|;i,0) 

(42) 

where P^{l2,X2\h,xi) and PAl2,X2\h,xi) are, respec- 
tively, the wedge and planar conditional probability of 
the interface being at a relative height I2 from the sub- 
strate at X2, provided that the interface is pinned at a 
relative height li at xi, defined as: 



Pi{l2,X2\ll,Xi 



Pi{h,xi;l2,X2) 
P^{h,xi) 



(43) 



where the subscript i indicates if this probability is con- 
sidered in the wedge (i = w) or in the planar (i = tt) 
geometry. 

In view of the identity between the wedge and pla- 
nar conditional probability distribution functions we first 
consider the of a planar substrate. The conditional prob- 
ability can be obtained as: 



P:il2,x\h,0) 



2 ''Z„{li,l2;x) 



(44) 



For contact binding potentials, Eq. (|44|l can be written 
explicitely as: 



Kih,x\h,o) = xl—c ^ 



+e' 



-212012 



Yj S(ii+i2-93:)'' 

e 2x 

2'KX 



-hS6'erfc 



2'KX 



{h +I2- Ox) 



(45) 



If li is very large compared with (1)^^ = 1/2E6', we 
can identify two different behaviors of P^(^2, 2:|Zi, 0) as 
a function of I2 (see Fig. If a; < h/O, the con- 

ditional probability is basically the free interface condi- 
tional probability that fluctuates around an average value 
(hix)) — h — Ox, with a standard deviation of the order 
of ^/xJY.. For X > h/O, the conditional probability be- 
comes the 1— point planar planar distribution function 
P-wih) = 2E6' exp(— 2E6'/2), completely uncorrelated to 
the value of li. The transition between the two regimes 
occur in an x interval around Xt = h/O which has a width 
of order of y^2li/Y,d^ = ^/xt^\. These results are con- 
firmed by the exact evaluation of the first moments of 
the conditional probability: 



{qnii,x) = / dh qp:{i2,x\ho) 



(46) 
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FIG. 3: Illustration of a typical interfacial configuration 
pinned at h ^ (Z)-^ for x = (thin continuous line). We 
have set E = 1 (it defines the length scale), 6 — 0.2 and 
h = 500. The thick continuous line corresponds to the con- 
ditional average profile {h)'^ {h , x) , and the dotted lines cor- 
respond to max(0, {l2}'^{h,x) ± 3fl {h,x)) where ^±{11, x) is 
the conditional roughness. Any interfacial configuration has 
a probability of at least 95% of being between the dotted 
lines. Inset: an enlargement of the area around Xt = li/6. 
Other characteristic length scales are represented. See text 
for explanation. 



The average conditional interfacial profile, which corre- 
sponds to n = 1, is given by: 



1 li-xO' 



I E _S(li~8x)f 
_ e 2:0 

ZTTX 



4E6I 



rfc 



4E6I 



-erfc 



2'KX 



{li+9x) 



(47) 



and the conditional roughness ^'j_{li,x) is defined as 
V m)'^ ~ {{h)'^)'^, where (Z|)'^(Zi,a;) can be written as: 



{iinh,x) 



ih - 9xf + 



^1 + 



{h - xOf 



erfc 



+ 



x9 -f h 



2E6I 



-erfc 




1 



2S 4E26/2 



2'KX 



ih - 9x) 



— {k+ex)\ (48) 



We obtain two main conclusions from these results when 
h ^ (Ott- First, the interfacial positions are highly corre- 
lated to the the central one for \x\ < li/O. Secondly, the 
intrinsic interfacial fluctuations are small in this x range 



compared to the conditional average value. Actually, if 
we set li as the length scale, the rescaled conditional 
probability distribution function P^{l2/li,x/li\l,Q) = 
hP^{hix\li,Q) behaves as: 

P-^{h/h,xlh\l,Q) ^ <5 f^l-^Ij±^\ H{h - ex) 



+5[f^j H{ex-h)m 

when YiOli oo. We expect this result to be valid for 
any potential and also for random bond disorder, since 
in all these cases the wandering exponent for the free 
interface C < 1- this can be checked for the marginal 
\/l potential. The Laplace transform of the conditional 
probability distribution as: 

/>oo 

C[P^[l2,x\hM = / dxe''^P:{l2,x\lu()) 
Jo 

= ^^z^{iij2,E-j:ey2) (50) 

For S — !■ oo at fixed E, 9, li and I2 , and taking into 
account Eq. lfTF|l and that the ground state eigenfunction 
ipoil) Of Z*^ exp(— S6'/), we obtain the following behavior 
for the Laplace transform of the conditional probability 
distribution function: 



1 



E(h~l2)/9 



5{l 



E 



C[P^{l2,x\h,0)] ^ -H{h-l2)e 

(51) 

The Laplace transform can be inverted, leading to Eq. 
(El. 

To proceed, we return to our discussion about the 
wedge geometry. Due to the presence of the boost fac- 
tor exp(2EaZ) in the midpoint probability distribution 
function, the midpoint interfacial height is almost always 
further from the substrate than the mean wetting layer 
thickness (1)^^ for any binding potential. If we assume 
that the conditional probability distribution function is 
given by Eq. H49|l , which corresponds to as neglecting the 
intrinsic interfacial fluctuations around the conditional 
interfacial proflle, we can capture the main features of 
both the average interfacial proflle and the correlations 
along the wedge for contact binding potentials. Actually, 
this picture is completely equivalent to the 2D wedge 
breather mode model 0, 0] . 

The average interfacial profile can be written as: 



{l{x)). 



dhPUh.O) 



dhPUh,O)ih~0x) 



dhhP^ih.^lluO) 



sP^{s + 9x,0)d452) 



The behavior of {l{x))u! for large x is dominated by the 
large I asymptotics of Pu,(Z,0). The latter can be ob- 
tained by taking into account Eq. H14|) for m = n = and 
making use of the WKB approximation for the 1— point 
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planar distribution function: 
1 / 



1 



2W(l) 



: exp 



e-2^'''exp(-2 / dt^] l^o^3) 



The first thing we can see is that, for large x, the de- 
cay of {l{x))u, in this approximation is controUed by an 
exponential term exp[—2Y,6{9 — a)x]. So, a new length 
scale is defined as 1/2'E9{9 — a). Close to the filling 
transition, = - IjTLaQ ^ (,f + C(l). 

Depending on the large I behavior of the (attractive) 
binding potentials, different situations can arise 0. The 
filling mean field regime is characterized by binding po- 
tentials that decay to zero as 1/1^ where p < 1/C — 1, 
implying C < 1 for thermal disorder (the wandering ex- 
ponent C = 1/2). A saddle point calculation shows that 
close to the filling transition {l{0))w ~ l/^(^ ~ ct)^- 
As — > a, the relevant length scale in the x direction, 
{l{0))w/0 » ^p, so the latter length scale is irrelevant (in 
fact, intrinsic interfacial fluctuations that we neglected 
can be more important). 

For p = 1, both length scales become of the same order, 
and consequently {l{x))io ~ {l{^))wf{x/£,p)exp{-x/£,*p), 
where f{x) diverges at most algebraically, and depends 
on the detailed structure of the binding potential through 
the short distance I dependence of Pw{l,0). For a pure 
l/l potential, f{x) = (1 -I- 2x/3 -|- x'^/6). This expression 
verifies the differential equation for {l{x))w that arises 
from Eq. H4l)|l in the scaling limit. 

The filling fluctuation regime corresponds to potentials 
with p > 1, and is characterized by universal critical 
exponents and scaling functions. Indeed in the critical 
regime the scaling behavior is the same as that found for 
contact binding potentials. For x ^ cxd, we find that 
asymptotically {l{x))u] ~ (Z(0))tu exp(— x/^|.). This solu- 
tion agrees with the asymptotics of {l{x))w for contact 
binding potentials when 9 ^ a, although with a decay 
length slightly smaller. However, the behaviour is asymp- 
totically correct if we assume that ^p = ^p. 

For the correlation functions, we have: 



{l{x)l{0))^~{l{x))^{l{0)),,= / dhhP^{h,0)A{h,x) 

Jo 

(54) 

where A(/i,x) is defined as: 

poo 

A{h,x)= dl2l2[P^{h,x\ll,0)-P^,{l2,x)] (55) 



In the breather mode approximation, A(/i, x) can be ob- 
tained as: 



A{h,x) « {h - 0x)H{h - ex) - {l{x))^ 



(56) 



We find different behaviors depending on the value of p. 
In the filling mean field regime, A(li,x) is negligible in 



this scale. For the filling fluctuation regime, the correla- 
tion function decays as: 



{i{x)i{o)),, - {i{x))^{i{o))^ ^ {imi 



1 



(57) 

and again is in agreement with the behavior of the exact 
correlation function for contact binding potentials Eq. 
(EH) when X — > oo and 9 a (assuming again that ^p = 
£,p). Finally, for the marginal case p — 1 the behavior 
of the correlation function is predicted to be for a; — s- oo 

{l{0))^g{x/S,p) eyip{—x/£^p), where g{x) is a function 
that diverges at most algebraically. 

Another quantity of interest is the midpoint local sus- 
ceptibility Xiu(Oj deflned as: 



Xi»(0 



dpii) 



dh 



h=0 



2{pi-p„) / dsP„(s,0)A(s)(58) 



where A{1) = dxA(l,x). In the breather mode ap- 
proximation and in the filling fluctuation regime, A(s) 
has the following expression: 



m)i 



(59) 



which is exact for contact binding potentials. This ex- 
pression, together with the midpoint wedge covariance 
Eq. |2Jl, leads to the covariance relationship between the 
local susceptibilities Q: 



-X^(^' 



(60) 



where d) is the local susceptibility corresponding to 
the planar geometry for a contact angle 9. 

Finally, we note that the breather mode picture has 
direct consequences for the scaling of the interfacial pro- 
file in the flUing fluctuation regime. To see this, recall 
that the wedge midpoint probability distribution func- 
tion scales as Q: 



P^ail) 



1 



-A 



I 



(61) 



where A(s) is a universal function and, due to covari- 
ance, is the same as the scaling function for the corre- 
sponding planar 1— point probability distribution func- 
tion. Complementing the scaling of the probability dis- 
tribution function is the position dependence of the in- 
terfacial profile, which we anticipate satisfies: 



{i{x))^ = nm^^cf) 



f Ox 



(62) 



where (j){s) is another universal function. In the breather 
mode picture, the interface is infinitely stiff in the filled 
region implying that the scaling functions A(s) and 0(s) 
are related via: 



<^(s) 



dsK{s){s - s) 



(63) 
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or equivalently: 



0"(s)=A(s) 



(64) 



A remarkable consequence of this relation is that the 
behavior of the interfacial profile close to the midpoint 
is determined by the short distance behavior of the 
wedge midpoint 1— point probability distribution func- 
tion. Since A(s) sc^^ as s — > we have (j){s) ~ 

1 — |s|-|-A|s| c for small s. Note that the first two terms are 
needed to preserve the continuity of the true interfacial 
profile {y(x))w and its derivative at the wedge midpoint. 
This result suggests that the interface behaves, for small 
values of x, as a random walk of a; as a function of z. This 
prediction is consistent with the behaviour of {l{x))w for 
contact binding potentials in the scaling limit 9 a, 
S.\\/S.F ^ but x/^p finite. 

Finally, note that Eq. (|64|l is also obeyed by the 
(non-universal) scaling functions corresponding to the 
marginal case. 



VI. RENORMALIZATION GROUP APPROACH 
TO THE CRITICAL FILLING TRANSITION 

In this Section we will justify the critical properties 
of the filling transition using a renormalization group 
framework. Specifically we will generalize an exact deci- 
mation functional renormalization group procedure pre- 
viously used to study 2D critical wetting 17, 18, 19] . 
Our transfer matrix results show that geometry plays a 
fundamental role in determining the critical behavior, so 
we anticipate that the appropiate renormalization group 
procedure must preserve the wedge shape. This implies 
that the effective wandering exponent ^ determining the 
rescaling of the interfacial height I must be C = 1- This 
contrasts with the value ^ = 1/2, which is appropiate 
for free interfaces and also planar wetting transtions. We 
will see that this choice leads naturally to the breather 
mode picture of the filling transition, implying that in- 
terfacial fluctuations are irrelevant except for those that 
determine the wedge midpoint interfacial position. 

Before introducing the renormalization group scheme, 
we generalize some of the results of previous Sections. 
The set of (2n -I- 1)— point distribution functions that in- 
cludes the midpoint interfacial position can be obtained 
in terms of the planar case counterpart by Eq. (|14|l . 
If we set 9 = a, it is clear from that expression that 
Pwi—n;...;n) = at the critical filling transition for 
any value of n. However, all the correlation functions 
decay at the same rate, since Eq. (|14|) can be rewritten 
as: 



P„(-n;...;n) = P„(0) 



P^(-n;...;n) 



(65) 



Consequently, the conditional (2n+l)— point probability 
distribution function remains finite at the filling transi- 
tion. The only relevant operator (in a renormalization 



group sense) should be related only to the 1— point prob- 
ability distribution function at the midpoint. Taking into 
account Eq. (|ll|l and the definition of the 2— point condi- 
tional probability distribution function Eq. 144f) . we can 
rewrite Eq. (|55|l as 



i=0 

II P^(k,x,+i~x,\k+i,()) (66) 



where we have chosen the ground state eigenfunction to 
be real and positive. The 2 —point conditional proba- 
bility distribution function has a non-trivial limit when 
Y,9li oo (see Eq. (gHJ)- Our goal will be to find 
a renormalization group scheme in which the 2— point 
conditional probability distribution function converges to 
this limit, and the only relevant operator is related to the 
wedge midpoint 1— point distribution function. 

Let us consider a discrete version of the interfacial 
hamiltonian Eq. |j2Jl: 



n-l 



PH = ^ I -(z,+i - z,)' + Wiz, - la) 



(67) 



where the spacing between sites a — 1 defines the length 
unit for I, S~^, etc. Using a similar transformation to 
the continuous case, Eq. H67|l can be written as: 



j3H[l] = 2n 



Ea2 



2'Saln — 2T,alo 



where periodic boundary conditions have been applied 
{In — l-n)- To simplify our discussion, we will consider 
n — > cx) and neglect the boundary effects. The probabil- 
ity of any interfacial configuration is given by: 



PM}) 



-/3H 



2T.alo+l3fw 



z 



n 



-PH{li,li+i) 



(69) 



where 7i.{liji+i) is defined as: 



■70) 



where fs is the planar surface free energy per unit length 
and is related to the contact angle corresponding to the 
binding potential W{1) via: 



2 



(71) 



For 9 = 0, Pfs converges towards the free interface free 
energy per unit length. Note that in the continuum limit 
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Eq. ©, this term is absorbed in the path measure. On 
the other hand, Pfw is defined as: 



/3/i 



w 



lim 



-In Z + 2n- 



(72) 



that corresponds to the wedge excess free energy (we 
suppose that this quantity is well-defined). Note that 
Ti.(li^li+i) is invariant under an exchange of its argu- 
ments. 

Let us now consider a decimation procedure, similar 
to the one used for the study of 2D critical wetting. We 
group the sites in blocks of b units, keeping the first one 
and integrating over all the interfacial positions of the 
remaining 6—1 sites in the block. Since the site i = 
plays a special role, the sites to be kept in each decima- 
tion step are those with i — jb, with j G Z. After that, 
we rescale x positions by a factor b and / by a factor 6^, 
i.e.: 



X 

b 



(73) 



The new hamiltonian PH'{ljJj^^) is defined as: 

poo 

Jo 



(74) 



Note that the renormalized hamiltonian /37t!'(/i, Z2) is 
symmetrical under an exchange of its arguments provided 
that the original /3Ti{li, I2) is also symmetrical (even if it 
is not defined as Eq. (jTOl). This procedure is iterated, 
leading to a sequence of renormalized hamiltonians. 

In order to complete the description of the renormal- 
ization group (RG) procedure we should give the trans- 
formation rules for a and f3fw- First we revisit the pla- 
nar geometry {a = fw = 0). We will consider the value 
of the exponent ( arbitrary, unlike in Refs. [T^ IT^. 
where ^ = 1/2. In general, after an arbitrary number of 
RG steps, we can write any j3'H{li, h+i) as: 



2 V27r/ 2 ^ ' 



where W{li^ U+i) is a symmetric function under exchange 
of its arguments. Obviously, in principle this function 
need not decay to zero when both arguments are large 
(as it does in Eq. ifTUI) '). However let us suppose that it 
decays as —A/[[li + lij^i)/2Y when kji+i — *■ 00. After 
making a RG step, we would like to find the asymp- 
totic behavior of the renormalized W'{l'j,lj_f_i) for large 
enough l'^, l'^^-^. Expanding the RHS of Eq. ((TH) and 
keeping terms up to first order in W (since I'j and I'jj^-^ 
are large, the values of Zi . . . Ib-i that contribute most to 



the integral are also very large), we find that: 



, rf^iV^^-P 2 



— exp 



27r " \ 2 / 

1-W{b^l'j,h)-W{k.i,b^l'^+,) 



b-2 



J2w{k,k+i) [76) 



i=l 



The lowest order term (corresponding to set W ~ 0) can 
be estimated for large I'j, h, . . . , h-i, I'j^i by extending 
the lower integration limits to — cx), and has the value: 



^ib<l'-hf 



fo'^e''^ / dhJ^exp 



dZ.-i^^exp 



e'-V 



2 /S52C-1 



27r 



■ exp 



Sfe2C-i(/'_/' 



So, in order that the renormalized binding potential W 
decays to zero at large values of both its arguments, the 
interfacial stiffness and the contact angle must transform 
as: 



S ^ E' = 



= eb^-^ 



(78) 



Two comments are pertinent at this point. Firstly, the 
transformation of the interfacial stiffness is also valid for 
a free interface. Secondly, the change in the contact an- 
gle has a geometrical interpretation, since its scaling Eq. 
(|78|l corresponds exactly to the change of small angles 
under the coordinates scaling Eq. (|73() . Thus it is rea- 
sonable to expect that a must change in the same way. 

If we take into account the first order in W , we can 
characterize the decay of the renormalized potential at 
large values of 1'^ and I'j+i- For simplicity, we consider 
the case 6=2. The renormalized binding potential has 
the following expression: 



— / dhe 



dlie 



Wiil'^,l'^+,) + Wi{l'^,l'^+,) (79) 



The first term Wl corresponds to the contribution of the 
hard wall to the renormalized binding potential, and W2 
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corresponds to the contribution of the original binding 
potential W. The hard wall contribution can be evalu- 
ated exactly as: 



1 



22CI] 



exp 



+ i-j+i 



(80) 



For W2, we take into account the long distance behav- 
ior of W. After some algebra, the leading order of 
W2{lj,lj^i) can be written as: 



w^2(n-,n+i) 



ds- 



{I'jT J -so/2 [l + s/soT 



(81) 



where l'^ = {V^ + Z;+i)/2 and sq = ^2^+^]. As 00, 
the integral tends to 1 and we can see that W ^ W2 ^ 
A'/{l'^)P, where A' = Ab^-'^P. It is interesting to note 
that this result is also obtained by the following scaling 
argument for the binding potential: 



A 



A' 



A 



Ab^~^P 



(82) 



{V)p "IP {bQ')P {l')P 

where we have taken into account that the binding po- 
tential is a free energy per x unit length. 

This analysis shows that the RG procedure leaves in- 
variant the functional dependence of the asymptotic be- 
havior of the binding potential. Two regimes can be 
identified. If p > 1/C, the binding potential strength 
decreases in each RG step. For p < 1/C the binding 
potential strength grows in each RG step. Finally, the 
marginal case p — 1/C, corresponds to the leading asymp- 
totic behavior remaining invariant. For C, — 1/2, which is 
the relevant value for planar wetting phenomena, this be- 
havior leads to the existence of two and three fluctuation 
regimes for complete and critical wetting, respectively 
|20l |. Note that this choice of C, keeps the relevant mi- 
croscopic scale ff, ^ invariant. The analysis of the 
critical wettiiig from this RG approach can be found in 
Refs. [l3,[ll|l|. 

A special class of effective hamiltonians are the follow- 
ing: 

exp[-/37Y(/„Z,+i)] = Zf^^k+i-^e-"^ (83) 

where Z'^'^[liji+i;l] is the partition function Eq. © 
with X — a = 1 for arbitrary values of the interfacial stiff- 
ness and binding potential W(l). Such Hamiltonians can 
be regarded as those which are generated after one iter- 
ation of the renormalization group provided that b very 
large. Indeed the fixed points found in Refs. p^lia.ll9l| 
belong to this class. Taking into account the properties 
of the path integrals, the renormalized potential after a 
RG step Eq. (|74ll can be written as: 



exp[ 



pn'{l',l'+,)]^b<Zf''^[b<l',b^l';b] 



(84) 



For contact binding potentials Eq. @ and Kratzer po- 
tentials Eq. H18|l . it can be checked that Eq. (|84|l corre- 
sponds to an effective hamiltonian of the same form as 
the original since: 



biz^-^[b%,b%+i;b] = zr''^'[l',,l',+,;i] 



(85) 



where S and 6 are transformed via Eq. (|78|l to S' and 
6'. For Kratzer potentials w must change as i.e. 
w' = wb^~'^'' = w/b in order to preserve the invariance of 
the leading order behavior under renormalization. 

Finally to finish our discussion of the RG for planar 
critical wetting phenomena, we note that the 1— point 
distribution function rcnormalizes as: 



(86) 



We will require this result later. Returning to our dis- 
cussion about the RG in the wedge geometry, we need to 
provide the transformation rules for a and Pfw ■ We will 
assume that a changes as 9: 



(87) 



In order to obtain the transformation rule for Pfw: we 
consider how the 1— point midpoint wedge probability 
distribution function rcnormalizes: 



p;(;^„o)^p;(;^,)e2^'"''o+(''/-)' 

= 6?P,(6?Z^,)e2^"''''o+(/9/™)' = b<PUb^l'o,0) 
implying that /3/^ remains invariant: 



(89) 



Finally, we note that if we change the effective hamilto- 
nian by: 

PH{k,k+i) = mik, k+i) ± [f{k+i) ~ fik)] (90) 

where the sign is positive for i > and negative for i < 0, 
the probability of an interfacial configuration is now: 



P^({?,}) = e2^"'°+2/('«)+''^"' n e" 



(91) 



The renormalization of the hamiltonian Eq. (|74f) is valid 
provided: 



f'{ll) = f{b^l'i) + C' 



(92) 



Note that any function /{h) (unless it is a constant) 
breaks the invariance of /37i under exchange of its ar- 
guments and consequently introduces a directionality in 
the X axis. This is perfectly sensible in the wedge geome- 
try, but is not admissible for the planar substrate, where 
—a; is completely equivalent to x. A convenient choice 
for f{li) is: 



1 



(93) 



13 



where PTr(li) is the 1— point probabihty distribution func- 
tion in the planar geometry, and the condition Eq. I|93|) 
is verified due to Eq. H86|l. ft is not difficult to see that 
exp(— /37i) is the 2— point conditional probability distri- 
bution function. Taking into account this fact and Eq. 
()88|l . it is clear that Eq. (|91|l is exactly the same as Eq. 
(PI- 

The one remaining issue to be decided is the relevant 
value of the exponent (. We choose C = 1, so that the 
wedge tilt angle a and contact angle 9 remain invariant 
in each step of the RG. 

The procedure is now standard. The RG flow trajecto- 
ries are constrained to the 6*— constant hypersurfaces in 
functional space. Since we know that the filling transi- 
tion occurs for 9 = a, we first check that this situation 
corresponds to the critical manifold. Instead of consid- 
ering an arbitrary potential, we choose as initial effec- 
tive hamiltonians those of the form given by Eq. (|83|l . 
in particular with the partition function corresponding 
to contact binding potentials and those of Kratzer form. 
When the number of RG steps n — > oo, the probability 
distribution of an interfacial configuration converges to a 
fixed point of the form: 



t3H'{li,h + i) 



(94) 



where {(ifw)* = — oo, and (30.* is defined as: 



{ 5{k+i -k + a)H{k 
-5{U)H{a-k) 



a) 



i > 



Sik- 
+S{h 



k+i + a)H{li+i - a) 
yi)H(a - k+i) 



i < 
(95) 

It is straightforward to see that Eq. (|95() is a fixed point 
of Eq. (|7H) . Obviously, since {(ifw)* — — oo, the proba- 
bility distribution for any interfacial configuration is zero. 
Nevertheless, we will consider Eq. (|94|) as formally dif- 
ferent from zero. The basin of attraction for the (j) — Q 
case is expected to be those hamiltonians with binding 
potentials that decay faster than \/l (corresponding to 
the filling fluctuation regime). Hamiltonians that have 
a binding potential with a leading order —(pa /I will be 
attracted to the fixed point Eq. (|94|l with the same value 
of (p for 9 = a (marginal case). Finally, if IW{1) diverges 
as Z ^ oo, there is no fixed point and the filling transition 
is mean-field-like. 

li 6 a, no matter how small \9 — a\ is, the RG flow 
drives the hamiltonian away from the critical manifold. 
For 6 > 0, exp(— /37i) converges to a fixed point ex- 
pression Eq. (|95|l with a replaced by 9. Even though 
this new fixed point for the conditional probability is dif- 
ferent from the the critical fixed point, the flow for the 
2— point conditional probability distribution function will 
remain close (in some functional sense) to the critical one 
if 16* — a| <C a. Consequently, there is no relevant field 
associated with PTi. 



On the other hand, the wedge midpoint 1— point prob- 
ability distribution function behaves differently. We con- 
sidered the same initial effective hamiltonians as for the 
6 — a case. If 9 < a, the distribution has an unphysical 
exponential growth with Iq as exp(S\9 — a\lo) . The expo- 
nential term grows in each RG step, driving the probabil- 
ity distribution function to infinity (however, Pf^ = — oo, 
so the "real" probability of any interfacial configuration 
is zero) . The attractor at infinity can be regarded as the 
complete filling fixed point. For 9 > a, the wedge mid- 
point 1— point probability distribution function becomes 
more and more peaked around zero as S — > oo, converg- 
ing to the low temperature fixed point: 



(96) 



These results imply that there is a relevant field (in the 
RG sense) associated with Ptu(/o,0). Actually, the only 
other relevant operator is /i oc fic — fJ-, where fi is the 
chemical potential and fj,c the value at gas-liquid coexis- 
tence. 

Recall that the critical exponents defined at coexis- 
tence (fi = lie): 



a(0))^ 



^^(0) 
I3fw- 



(97) 



where t = Tf—T and Tf is the critical filling temperature. 
Close to a critical filling fixed point (for any (/>), all the 
relevant scale lengths (^(O))^,, £,x, £,±{0), etc. are 

reduced by a factor 1/b. To extract the dependence on t 
we need to know the largest eigenvalue of the linearized 
RG flow close to the critical fixed point. We argued above 
that this eigenvalue is associated to the transformation 
of the a; = term of the hamiltonian PHq. We again 
considered the special hamiltonians Eq. (|83|) for contact 
and Kratzer binding potentials. The expression of (3Hq 
for these potentials is: 



PHq = 2E(0 - a)lo - 20 In ?o + C 



(98) 



Taking into account how f3Ti.Q renormalizes, and its ex- 
pression at the critical fixed point (see Eq. H94() '). it is 
clear that the largest eigenvalue is b, and its associated 
eigenfunction is proportional to E(0 — a)l(). As a con- 
sequence, the critical exponents i'± = Vx = Pw = ^ for 
both the filling fluctuation and marginal regimes. 

To get the leading singularity for the wedge excess free 
energy is more complicated. As (if^ remains invariant 
under the RG transformations, 2 — = 0. However, 
this value does not rule out a logarithmic divergence in 
— T. In order to get such a dependence, we use the 
following thermodynamic relationship: 



d{f3U) 
da 



= -2E(?(0)), 



(99) 



s.e 



where the derivative is made without changing E or any 
characteristic of the binding potential (in particular the 
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contact angle). We can rewrite Eq. (|99|l as: 

After a renormalization step, we know that, even when E 
and the binding potential have changed, this derivative 
changes as \/h since /3/u,, 9 and a are invariant, and 
and decrease by a factor 1/h. If we regard 

Ti{0 — a) (oc S, and recall that 9 — a is invariant under 
a RG step) as proportional to t when the flow is close 
to the critical fixed point, Eq. H1UU|I implies that (3f^ 
diverges logarithmically as T Tf . Furthermore, Pfw 
is invariant under renormalization and must vanish as 
a — > 0. Consequently close to the filling transition it 
must be proportional dependence to hi(9 — a) — hi(9). 

The critical exponents obtained for the critical filling 
transition are in complete agreement with exact calcula- 
tions and scaling arguments (^]. 



breather mode picture of fluctuation effects. The same 
picture also emerges of a renormalization group approach 
that leaves the wedge geometry invariant. The fact that 
the only relevant fluctuations are those that translate the 
midpoint interfacial height leads to a simple relationship 
between the interfacial structure and midpoint height 
probability distribution function in the scaling limit. The 
extension of this approach to higher dimensions and/or 
different geometries would be very interesting and further 
work is being carried out in that direction. 

Finally it would be instructive to understand the phys- 
ical origin of the covariance relationship Eq. (|2Jl. In the 
filling fluctuation regime we found that covariance can be 
inferred from the existence of the differential relation Eq. 
(|25ll and some regularity conditions. The very existence 
of such a field equation is itself indicative that some un- 
known symmetry relates wetting and filling transitions. 
Further work is required to elucidate whether any such 
hidden symmetry exists. 



VII. CONCLUSIONS 

The structure of the gas-liquid interface bound at a 2D 
wedge and close to the filling transition has been stud- 
ied using exact transfer-matrix methods. These calcula- 
tions show the emergence of a new correlation length 
sufficiently close to the phase boundary {a < 9 < 2a). 
The explicit transfer matrix results for correlation func- 
tions and interfacial roughness completely support the 
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